Near-field radiative heat transfer between macroscopic planar surfaces 
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Near-field radiative heat transfer allows heat to propagate across a small vacuum gap in quantities 
that are several orders of magnitude greater then the heat transfer by far-field, blackbody radia- 
tion. Although heat transfer via near-field effects has been discussed for many years, experimental 
verification of this theory has been very limited. We have measured the heat transfer between two 
macroscopic sapphire plates, finding an increase in agreement with expectations from theory. These 
experiments, conducted near 300 K, have measured the heat transfer as a function of separation 
over mm to /im and as a function of temperature differences between 2.5 and 30 K. The experiments 
demonstrate that evanescence can be put to work to transfer heat from an object without actually 
touching it. 

PACS numbers: 44.40. +a,78. 20. Ci 
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Humans knew of radiative heat transfer at least as 
early as the discovery of fire, and physicists have investi- 
gated this process for centuries, culminating in the black- 
body theory of Planck and the birth of the quantum the- 
ory. Planck's equation for black-body radiation contains 
only the temperature and some fundamental constants. 
When actual materials are involved, their emissivities en- 
ter the discussion, but little else. For example, the heat 
transfer per unit area between two semi-infinite planes 
is set by their temperatures and integrated emissivities 
but does not depend on their separation or other geo- 
metrical quantities. When the two planes approach each 
other closely the situation changes. In this near-field 
regime, each material interacts with exponentially de- 
caying evanescent electromagnetic fields generated in and 
existing outside the other material; these fields can drive 
currents and generate heat.[lr|3| This near- field radiative 
heat transfer can be several orders of magnitude greater 
than far-field blackbody radiation. 

Much like the Casimir and van der Waals force, near- 
field heat transfer deals with fluctuations that only exist 
over small distances. The first in-depth theory for near- 
field heat transfer between planar surfaces was derived by 
Polder and Van Hove,@] building on the work of Rytov[lj . 
There have been several other theoretical approaches, 
and in general the theory seems complete, except per- 
haps at distances comparable to atomic dimensions. [4{ 

Although heat transfer via near-field effects has been 
discussed for many years, experimental verification of the 
theory for heat transfer between two planar surfaces has 
been limited. HargreavesfBj has presented room tem- 
perature observations for two Cr surfaces at distances 
as small as 1 /im. Domoto et al.@ reported results at 
cryogenic temperatures but for relatively large (50 /xm) 
separations, where near-field effects were barely observ- 
able. Neither study compared experiment to theory. A 
comparison at a fixed spacing has been put forward, but 



the plates were separated by polyethylene spacers, so the 
distance could not be varied. [7| There have also been sev- 
eral recent results using a sphere-plane geometry. 0, 
There are engineering reasons for this approach, as, un- 
like a parallel plane geometry, a sphere-plane geometry 
needs no angular alignment. Using scanning probe and 
micro-machine technologies, these experiments cover a 
wide distance range. 

In this paper, we report a set of measurements of heat 
transfer between two planar surfaces and present a de- 
tailed comparison to theory. In addition to its intrin- 
sic interest, this work was motivated by possible ap- 
plications in cooling objects without actually touching 
them, such as the mirrors of a future laser interferometer 
gravitational- wave detector, [lfj This application would 
require the parallel plane geometry in order to get large 
areas of close approach of the two objects and thus sig- 
nificant heat transfer. 

The near-field heat-transfer process can be thought 
of as a form of frustrated total internal reflection. 
Evanescent waves, exponentially-decaying electromag- 
netic fields, exist outside, but near to, the surface of 
a material medium at temperature T. These decaying 
fields are a consequence of travelling waves inside the 
medium experiencing total internal reflection; this phe- 
nomenon occurs when there is no valid solution to Snell's 
Law, rii sin Qi = n t sin 9 t , where rii and rit are the indices 
of refraction and 6$ and t the angle between wave vec- 
tor and surface normal on the two sides of the interface. 
Although there is no energy transmission across the inter- 
face, an electric field exists on the far side. Furthermore, 
if another medium is brought near this exponentially de- 
caying field some of the energy from the incident beam 
will propagate across the gap and into the new mate- 
rial. If one medium is hotter than the other, this photon 
tunneling will lead to heat transfer from hot to cold. 

Polder and Van HoveQ considered two half spaces of 
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identical material but different temperatures. The com- 
puted heat transfer coefficient W comes from the tem- 
perature derivative of the z component of the Poynting 
vector (S z ) from each medium, which is evaluated just in- 
side surface of the other medium. W is the temperature 
derivative of this difference as can be seen in Eq. [TJ 

S z (d+) - S z (0-) dS z (d+) 
Ti-T^o 1 Ti-T 2 1 dTi K ' 

It is convenient to break W into two components 

W = W sin + W cxp . (2) 

Wsin is the part where the wave number k x of the field 
is in the range < k x < lu/c, where the field is propa- 
gating, and which gives the ordinary (Stefan-Boltzmann) 
far-field radiation. W cxp is the component where the field 
is exponentially decaying away from the surface; it is the 
larger contributor to the near-field limit: 
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where the energy transmission coefficients are 

1 - cos(2xp) 
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cosh[2/«(<i — dp)] — cos(2xp) ' 
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with P = || or _L for the two polarizations and \\\ an d 
X± being the phase shifts on reflection 

X|| = arg[(-i«e + k z )(-ine - k g )*], (5) 



X± = axg[(-iK + k z )(-iK-k z )% (6) 

with e the material dielectric function. 8\\ and 5j_ can be 
derived from 
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and k z and k are the z component of k for the medium 
and the vacuum, respectively 



k z = y/(e - l)^ 2 /c 2 



k = ik zv = y/k% — U) 2 /c 2 . 



(9) 



(10) 



Note that k x > w/c, so that k zv is imaginary above the 
surface and the energy transmission coefficients 7jj° xp and 
7^ xp arc not derived from the energy reflection coeffi- 
cients 1Z\\ and 7?a.[2| We have used this theory to calcu- 
late the heat transfer coefficient between sapphire plates, 
using Barker's sapphire dielectric function. 




Separation (fim) 

FIG. 1. Heat transfer coefficient vs. distance for z-cut sap- 
phire. The temperatures of the two media are T^ot = 310 
K and T C oid = 300 K, respectively. The violet curve shows 
W]- n + ITsin which dominate at far distances. The red curve 

shows wJxp. The green curve shows W^ p . The evanescent 
terms dominate at close distances. The blue curve is the total 
heat transfer coefficient W. 



Figure [T] shows the prediction of the model for z- 
cut sapphire; the temperatures of the two media are 
T hot = 310 K and T cold = 300 K. The W sin term dom- 
inates at large separation and is nearly constant with 
separation. The W eX p term dominates in the near-field 
regime, so that the total heat transfer changes from be- 
ing independent of separation at large distances towards 
W cxp oc l/d 2 at short distances. The turning point be- 
tween Wsin and Wcxp occurs when the separation is ap- 
proximately equal to the peak wavelength of the black- 
body curve of the hot half space. Wien's displacement 
law predicts that X ma x ~ 9 fim when the temperature is 
310 K. 

A sketch of our apparatus is shown in Fig. O It is de- 
signed around two 50x50x5 mm 3 sapphire plates. These 
have a specified flatness of A/8 @ 633 nm per inch on the 
largest surfaces and are cut such that the c axis is per- 
pendicular to these surfaces (z-cut). Sapphire was used 
because it has good thermal conductivity. It is also a 
candidate for the test masses of future gravitational- wave 
detectors. 

One of the plates, henceforth called the cold plate, is 
attached to the thermal bath that is the vacuum cham- 
ber. The other plate, the hot plate, is thermally isolated 
from the bath by a Macor spacer attached to the back 
side of the hot plate. The hot plate also has a heater 
wound on a copper ring which is itself attached to the 
back of the plate. The heater current and voltage, after 
correcting for lead resistance, give the power required to 
maintain a given temperature difference between the hot 
and cold plates. Both plates have a Si-diode thermome- 
ter fastened to their backs to read the temperature (and 
for the hot plate to control it). Both plates have all four 
corners coated with an approximately 200 nm thick layer 
of sputtered copper. These coatings have areas about 



3 



Mirror Mount for 
tip/tilt/translation 
Heater control ~~] 



■motor (x3) ,„ nfrnMi 



motor 
controller{x3) 




Thermometer 



FIG. 2. Experimental apparatus. Stepper motors allow ad- 
justment of the spacing, tip, and tilt (read capacitively) of two 
sapphire plates. The temperature of the hot plate is controlled 
by a feedback circuit, and the power required to maintain a 
temperature difference gives the heat transfer from the hot 
plate to the cold plate and to the thermal bath. 



1 mm 2 and serve as capacitor plates that are read by 
four 24-bit capacitance-to-digital converter circuits [12| to 
measure the separation and angular misalignments of the 
plates. The metal film is wrapped around to the sides of 
the sapphire to allow electrical contact to the electrodes. 

The cold plate is glued to a copper disk, which in turn 
is attached to the experimental structure. The Macor 
spacer on the back of the hot plate is attached to a modi- 
fied kinematic mirror mount which allows for z-axis linear 
movement and tip and tilt angular adjustment by turn- 
ing the three adjustment screws in the back. Three step- 
per motors turn screws on the kinematic mount via gear 
reduction boxes; each motor step translates to a linear 
movement of the hot plate by 35 nm. The components 
are held together by an "L" shaped backbone (not shown) 
to give rigidity. The assembly is located in a UHV cham- 
ber, with a base pressure below 2 x 10~ 7 Torr, making 
gas conduction negligible. Signals to the stepper motors, 
capacitance readouts, temperature readouts, and current 
and voltages to the heater are all controlled and/or read 
by a Lab VIEW computer program. 

Each pair of capacitor plates is calibrated by taking ca- 
pacitance readings as the plates are driven together one 
step of the stepper motors at a time. A fit is made to 
C = eoa/d+ C s tray where eo is the dielectric constant of 
the vacuum, a is the capacitor area, and C s t r ay is a par- 
allel contribution independent of separation. The data 
fits the equation above very well, with R 2 values greater 
than 0.999. The fitted value of a equals the metalized 
area within our knowledge of this area; C s t ra y ~ 0.4 pF. 
The average capacitance gives the distance while the in- 
dividual readings are used to correct the alignment by 
sending steps to the motors controlling tip and tilt. 

To measure the heat transfer coefficient we compute 
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FIG. 3. Heat transfer coefficient vs. distance. The curves 
are each offset vertically by 2 W/m 2 -K from the one below. 
The points are the data, with error bars, determined from the 
scatter in the heat transfer measurements and the uncertainty 
in the distance calibration. The solid lines are the theoretical 
predictions for flat plates while the dashed lines are the theo- 
retical predictions for slightly curved plates (see text). Each 
measured curve has a reproducible addendum due to other 
heat leaks which are not included in the model and which has 
been subtracted from the data. The temperatures are (top to 
bottom): T hot = 327.0 K, T co i d = 308.0 K; T hot = 322.0 
K, T cold = 307.0 K; T hot = 317.0 K, T co i d = 305.8 K; 
T hot = 312.0 K, T cold = 305.2 K. 



W = P/[A(T hot - Tcoid)], where P is the power dissi- 
pated in the heater, A is the plate surface area, and Thot 
and T C oid are the temperatures of the hot and cold plates, 
respectively. The data are a sum of parallel heat path- 
ways, including thermal conduction through the Macor 
spacer and the other parts of the hot-plate holder, radi- 
ation to the thermal bath, and the contributions of the 
near- and far-field radiation between the two plates. We 
can observe the near-field effect because it is the only 
one of these that will change with plate separation. The 
thermal conduction and radiation-to-the-bath paths just 
add a constant offset to the radiative heat transfer that 
the model predicts. 

Data for the heat transfer coefficient versus distance 
were collected for four temperature differences. (Sec 
Fig. El) Each one shows that near- field heat transfer ex- 
ists and that the data and model agree reasonably well. 
Each run covers a separation range of about 2-100 fim. 
The only freedom in the fit is an offset to the model. All 
our measurements have an offset of 0.0435±0.0004 W/K, 
completely consistent with thermal conduction through 
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FIG. 4. Near-field heat flux vs. distance for multiple temper- 
ature differences and multiple runs. (One run is shown with 
open symbols and other runs with filled symbols.) The inset 
shows the near-field heat flux vs. temperature difference for 
several specific separation distances. 



the Macor spacer and the rest of the support for the hot 
plate, in parallel with radiation from the rear surface of 
the plate. The chamber pressure for these runs ranged 
between 5 x 10~ 7 Torr and 2 x 1CT 7 Torr. The vac- 
uum chamber is held at constant temperature of 30.0 °C 
(303.2 K). The hot plate is brought close to the cold 
plate step by step. Each datum is an average of a set of 
500 values, in both heat transfer coefficient and distance, 
taken after the system has reached thermal equilibrium. 
It takes 30 to 45 minutes after moving to a new position 
to reach thermal equilibrium. We calculate the average 
and its standard deviation and plot these in Fig. [3l 

Theoretical curves match the experimental data well. 
However, although the agreement is within experimen- 
tal errors, there does appear to be a systematic offset: 
the theory predicts a slightly lower heat transfer coeffi- 
cient at each separation than we measure. Alternatively, 
the plates could be slightly closer than measured by the 
capacitive readout. We believe that the latter explana- 
tion is correct. Simulations of the heat transfer between 
two convex plates (shown as dashed lines) eliminate the 
systematic error when the radius of curvature in each 
plate is ~1 km, corresponding to deviations from flatness 
of 500 nm. Subsequently, we measured the plate cur- 
vatures using optical fiats and green mercury light; the 
Newton's rings interference pattern gave a central dis- 
placement of 170 ± 30 nm with respect to the perimeter, 
a little smaller but of the same order of magnitude as the 
simulation gave. 

Figure @] shows the heat flux caused by near-field ef- 
fects. Both far-field heat transfer and the offset due to 
heat leaks to the thermal bath have been subtracted. 



For each temperature difference, data are shown for sev- 
eral distinct data runs; these agree very well. The inset 
shows the dependence on the temperature difference of 
the plates. For distance values where the near-field ef- 
fects dominate, the heat flux is linear in the temperature 
difference. Both data and model were fitted to 
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where (p(AT, d) is the total near-field heat flux, G(d) 
is a multiplicative factor, and a(d) is the exponent for 
AT. Each curve follows a power law in AT, with the 
exponent varying from 0.70 at small distances to 0.91 at 
larger distances. The differential heat transfer, Eq.Q] has 
a(d) = 1; the finite temperature differences used in the 
experiment bring in higher-order terms. 

In summary we have measured near-field heat trans- 
fer across a small gap for a parallel-plane geometry. The 
data agree quite well with the theory of Polder and Van 
Hove[2j. The experiments demonstrate that significant 
amounts of heat can be transferred via evanescent radi- 
ation in the near-field regime. 
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